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A class of non-local contact processes is introduced and studied using mean-field 
approximation and numerical simulations. In these processes particles are created 
at a rate which decays algebraically with the distance from the nearest particle. 
It is found that the transition into the absorbing state is continuous and is char- 
acterized by continuously varying critical exponents. This model differs from the 
previously studied non-local directed percolation model, where particles are created 
by unrestricted Levy flights. It is motivated by recent studies of non-equilibrium 
wetting indicating that this type of non-local processes play a role in the unbinding 
transition. Other non-local processes which have been suggested to exist within the 
context of wetting are considered as well. 



2 



I. INTRODUCTION 



The contact process (CP) is known as a simple model for epidemic spreading that mim- 
ics the interplay of local infections and recovery of individuals It is defined on a 
(i-dimensional lattice whose sites could be either active (infected) or inactive (non-infected), 
denoted as '1' and '0', respectively. The model evolves random-sequentially by two compet- 
ing processes, namely, nearest-neighbor infections 01/10 11, and spontaneous recovery 
1^0. Depending on the relative frequency of these moves the contact process displays a 
continuous phase transition from a fluctuating active state into an absorbing state which 
belongs to the universality class of directed percolation (DP) 

In the present paper we investigate a generalized version of the one-dimensional contact 
process where inactive sites can be activated over long distances. We consider a lattice model 
where sites could be either active or non-active. The model evolves by random sequential 
updating with the following transition rates: 

1 with rate 1 (1) 
1 with rate q/t (2) 

The first move corresponds to the usual local annihilation process, while the second move 
describes a long range process in which an inactive site which is located at a distance / from 
the nearest active one becomes active. The rate at which such a process takes place decays 
algebraically with the distance I from the nearest active site (measured in lattice site units), 
reflecting the long range nature of the interaction. The overall rate of this process is governed 
by the control parameter q, while the characteristic shape of the interaction is controlled by 
the exponent a. In particular, usual short range dynamics is recovered in the limit a ^ oo. 
In what follows we refer to this model as the a-process. Note that its definition is valid in 
any dimension although we are primarily interested in the one-dimensional case. 



The long-range contact process investigated here is motivated 




jy recent studies of depin- 
lll. where one considers a 



ning transitions in non-equilibrium wetting processes 
fluctuating interface next to a hard-core wall. Regarding the pinned domains of the inter- 
face as active sites and unpinned domains as inactive ones, the dynamics of the fluctuating 
interface may be projected onto that of a contact process (see Fig. P). For example, unpin- 
ning of the interface due to deposition would correspond to spontaneous recovery 1^0, 
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FIG. 1: Long-range infections in the generalized one-dimensional contact process. Any site between 
two active sites can be activated with a rate proportional to where / is the distance to the 
nearest active site. The model is motivated by the dynamics of a fluctuating interface growing on 
an inert substrate (see text). 

while lateral growth of a pinned domain would correspond to local spreading of activity 
by 01/10 ^11. However, as unpinned regions of the interface fluctuate in the bulk they 
may spontaneously become pinned to the substrate far away from other binding sites. As 
illustrated in the figure, such a pinning effectively leads to long-range infections in the cor- 
responding contact process. Clearly, the rate for pinning will decrease as we move away 
from the binding site. Numerical studies of some wetting models (not shown in this work) 
indicate that this rate decays algebraically as with an exponent a ~ 2.6. This motivates 
the postulated power law in Eq. Nevertheless the relation between non-equilibrium wet- 
ting and the a-process should be considered as an analogy, rather than an exact mapping, 
because the projection neglects the internal dynamics of unpinned regions of the interface. 

In recent years the contact process has also been generalized in various other ways to 
include long-range infections Ij^Il^- In these studies it is assumed that a given active site 
can activate any other site with a probability P{k), which decays algebraically with the 
distance k as 

P{k) ~ k-". (3) 
In such models the infection can be thought of as being spread by Levy flights Q| . Janssen 



et al. solved this problem using field-theoretic renormalization group 



results of this analysis is in agreement with numerical simulations 
studies have been extended to contact processes with temporal 
Levy flights jlS,] as well as to systems with disorder |19| . 



calculations [l^. The 
1£ I . More recently these 
17( 1 and spatio-temporal 
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The a-process studied in the present work differs from Levy-ffight-mediated contact pro- 
cess. Here an inactive site can only be activated by the nearest active site. Therefore, in 
one space dimension, the interaction range is effectively cut off by half the actual size of 
the corresponding island of inactive sites. Thus infections in the a-process are mediated by 
truncated Levy flights, which cannot overtake other active sites. As will be shown below, 
this truncation changes the critical properties at the transition significantly. In particular, 
although both models exhibit critical exponents which vary with a, the exponents them- 
selves are not the same. Moreover, while in the Levy flights model the exponents become 

n 

mean-field like for small a jj] no such regime exists in the a-process. 

In the following section we analyze the a-process by a mean-field approximation as well 
as numerical simulations and compare it to contact processes with unrestricted Levy flights. 
In Sect, mil we further generalize the model by taking the rate for nearest-neighbor infection 
to depend on the distance to the nearest active site. This rate is controlled by an additional 
exponent a. The phase diagram of this generalized model, called a-cr-process, is then studied 
by mean-field and numerical methods. Conclusions are finally presented in Sect. IIVI 



We now consider the a-process within a mean-field approximation in terms of the density 
of active sites p{t). To this end the mean-field equation of the ordinary contact process 
has to be extended by a term that accounts for long-range infections according to Eq. 
Assuming that the sites of a one-dimensional lattice are uncorrelated and independently 
active with probability p, it is easy to check that the probability of an inactive site to lie 
at a distance I from the nearest active site is (1 — (1 — p)^)(l — p)^'~^- Summing up the 
contributions for all distances /, the mean- field equation governing the dynamics of the 
density of active sites takes the form 



II. ANALYSIS OF THE MODEL 



A. Mean-field approximation 




(4) 



1=1 
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where the first term appearing on the r.h.s. of Eq. ^ corresponds to the short range 
annihilation process. Note that the usual DP short range activation term is given by the 
first term of the sum in the r.h.s. of Eq. (j3)) 



Turning the sum into an integral the leading terms in p in the dynamical equation (ji]) 
may be evaluated. It is found that for a < 1 and, to leading order in p, Eq. (jlj) becomes 

dtp = rp + up"' (5) 

where (r, u) are constants and u > 0. Since the leading term in this equation is and since 
its coefficient is positive, the absorbing state (p = 0) is always unstable and no transition 
takes place for any finite q. This is also the case for a = 1, where logarithmic corrections to 
the leading linear terms destabilize the absorbing phase. 

On the other hand, for 1 < a < 2 the leading term in the equation is the linear one. 
Moreover the coefficient u of the leading non-linear term is negative. This results in a 
continuous transition to the absorbing state which takes place at r = corresponding to a 
non-vanishing rate q. Such a phase transition belongs to a universality class different from 
the one of the short range DP model. For example, in the stationary state the density of 
active sites scales as 

P.. ~ rV("~i) , (6) 

hence the mean- field order parameter exponent, = varies continuously with a. It 
diverges in the limit a ^ 1. For a > 2, however, the leading non-linear term is vp'^ where v 
is a negative constant. The mean- field equation thus reduces to that of the short range DP 
process, and usual DP-like transition is expected. 

Thus 

f ^ 1 < a < 2 
1^ 1 a>2. 

Before discussing the results of the numerical study of this model we consider a more 
refined mean-field approximation, where the effect of the long-range process on the diffusion 
constant is taken into account. To proceed we assume that, although the dynamical process 
is long range, it does not modify the usual Laplacian form of the spatial interactions, DV^p. 
The reason is that the infection range in the a-model is cut-off by the maximal distance to 
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the nearest active site, making it effectively of finite range. However, as we shall see below, 
the effect of the a infection process is to make the diffusion coefficient D dependent on the 
density p. 

In order to evaluate -D(p) we note that the infection process is cut off at length-scales 
^av ~ 1/ P, which is the average distance between active particles in a system with density p. 
The mean square displacement of an infection is expected to scale as 

~ / ~p°-3 + Do. (8) 

where A is a cutoff due to lattice spacing from which the usual constant contribution Dq 
(positive for a > 3) to the diffusion coefficient arises. Since -D(p) ~ (/^), Eq. (jH)) suggests 
that the effective diffusion constant diverges for a < 3. Thus the space-time-dependent 
version of the mean-field equation for 1 < a < 2 reads to lowest order 

aip(x, t) = rp{^, t) + n[p(x, t)r + [p(x, t)]"-3vV(x, t). (9) 

For this equation, dimensional analysis yields the complete set of mean-field exponents 

MF _ ^MF _ _J_ ^ = 1, (10) 

a — 1 " 

where as usual and are the mean-field values for the spatial and temporal exponents 
which control the power-law divergence of correlations (respectively in space and time) at 
criticality. 

On the other hand for 2 < a < 3 the non-linear term in the dynamical equation wp" 
becomes vp"^ while the diffusion coefficient remains singular at p = 0. This implies that the 
DP regime sets up only for a > 3, and not for a > 2 as suggested by the analysis of the 
reaction terms alone. In particular, for 2 < a < 3 lowest order mean-field reads 

dtp{^, t) = rp{^, t) + i;[p(x, t)]2 + [p(x, t)]°-3v2p(x, t). (11) 

and by dimensional analysis one gets 

^MF _ 2 - I , = z/^^ = 1 . (12) 

For the dynamical exponent z = we obtain 

a -1 l<a;<2 
= < 4^ 2 < a < 3 (13) 
2 a > 3 

In the following we compare these values with the results of extensive numerical simulations. 
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B. Numerical results 



In order to determine the critical exponents numerically, we performed Monte- Carlo sim- 
ulations of the one-dimensional a-process with periodic boundary conditions, starting with 
a fully occupied lattice. The time step dt of a single move in the random sequential updating 
procedure has been taken as 1/L, where L is the lattice size. 

For models with long-range interactions such simulations are difficult to perform because 
of strong finite-size effects, which increase as a — 1. In the present model these finite-size 
effects manifest themselves by sudden transitions into the absorbing state. This leads to 
uncontrollable fluctuations in the statistical averages. To circumvent this problem we com- 
bined flnite size scaling analysis with numerical simulations. The simulations were carried 
out on large systems with homogeneous initial states approaching the critical point from the 
active phase. 

The critical point q^. and the dynamical exponent z have been determined by flnite size 
analysis of the average absorbing time r for an initially fully occupied lattice. At the critical 
point r is expected to scale with the system size L as 

T ~ L\ (14) 

For one dimensional systems belonging to the DP universality class the dynamical exponent 
is known to be = 1.580745(10) j2l|. We also measured the exponent 5 = which 
is associated with the temporal decay of the density p of active sites at criticality. This 
calculation has been performed by applying the method of flnite-size analysis proposed 



by de Oliveira and Dickman [2^: whenever the system falls into the absorbing state the 
dynamical process is recovered by automatically resetting the system into a "typical" active 
conflguration. In practice, this procedure suppresses the large fluctuations which result from 
the abrupt transitions to the absorbing state in flnite systems. A simple scaling argument 
shows that the stationary density at criticality should scale with the system size as 

Pstat ~ L-''. (15) 

Finite size analysis has been performed by averaging over 10^-10^ independent runs. 

For a > 2.0 the estimates of the exponent 6, obtained by flnite size analysis, have 
been compared with the standard analysis, based on the time-decay of p{t) at criticality. 
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FIG. 2: Results of the numerical analysis of the a-process for a = 1.5. All plots are 
in double-logarithmic scale. a) Density of active sites p{t) as a function of time mea- 
sured for a system of size L = 2"^^ starting from a fully occupied lattice. In the ac- 
tive phase and close to the critical value Qc = 0.4695(5) (from top to bottom q = 
0.679, 0.569, 0.519, 0.489, 0.477, 0.473, 0.472, 0.471, 0.470) finite-size fluctuations yield sud- 
den transitions to the absorbing state, b) Finite-size scaling analysis of the average absorbing time 
r. From top to bottom q = 0.475, 0.471, 0.470, 0.469, 0.468, 0.467, 0.45. c) Finite-size scaling 
analysis of the stationary density pstat (see text), d) Stationary density pstat as a function of the 
distance from qc ■ The dashed lines in panels c-d are the best fit (see tabled). 



Usually, large fluctuations make the latter method quite inaccurate and time consuming. In 
order to obtain better performances we carried out numerical simulations of large systems 
(L = 2^'^ ~ 2^^) and over long time lapses (10^ ~ 10^ time steps), while averaging over a few 
different realizations. Good agreement with the predictions of finite-size analysis is found. 

The critical exponent (3 has been determined by measuring the stationary density of 
active sites p ~ {q — Qc)^ in large systems (L = 2^^ ~ 2^^) for different values of q, just above 
qc- 

The estimated values for the exponents 6 and (3 can be compared with the best numerical 
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b 


/3 






1.2 


0.205(3) 


0.34(2) 


2.0(1) 


2.73(15) 


5 


0.2 


1.5 


0.4695(5) 


0.67(2) 


0.84(4) 


1.25(5) 


2 


0.5 


1.8 


0.714(1) 


0.99(5) 


0.46(3) 


0.68(7) 


1.25 


0.8 


2.0 


0.8592(2) 


1.21(2) 


0.31(1) 


0.49(2) 


1 


1 


2.1 


0.9250(3) 


1.25(2) 


0.28(1) 


0.46(2) 


1 


1.0526... 


2.3 


1.0453(2) 


1.43(2) 


0.22(1) 


0.36(1) 


1 


1.1764... 


2.5 


1.1492(2) 


1.48(4) 


0.19(1) 


0.32(1) 


1 


4/3 


2.6 


1.1955(2) 


1.54(3) 


0.175(10) 


0.30(1) 


1 


1.4285... 


2.7 


1.2381(2) 


1.56(2) 


0.17(1) 


0.294(8) 


1 


1.5384... 


3.0 


1.3470(2) 


1.58(2) 


0.166(8) 


0.278(8) 


1 
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TABLE I: Estimates of the critical points and exponents for various values of ot. These values are 
plotted in Fig. 1^1 For comparisons mean field values are also reported. 

estimates of the DP exponents in one dimension, 5°"" = 0.159464(6) and (5°'^ = 0.276486(8) 
[2^. Typical curves resulting from this numerical procedure are shown in Fig. |21for a = 1.5. 



As summarized in Table HI the numerical results show a qualitative agreement with the 
mean-field predictions. In particular, it is found that all exponents vary as a is increased 
up to a critical value a, beyond which they become consistent (within numerical error) with 
the DP values. As a approaches 1, the exponents (3 and 6 are found to increase, while the 
dynamical exponent z vanishes as a — 1. The comparison between numerics and mean- field 
is shown in Fig. El 

Note that according to numerical results the crossover to the DP scaling regime seems to 
take place at a ~ 2.7. On the other hand, the mean-field predicts that the crossover takes 
place in two steps: first, (3 and 6 become DP-like at a = 2 and then the exponent z assumes 
its DP value at a = 3 . 
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FIG. 3: Comparison between numerical estimates (including error bars) and the mean-field pre- 
dictions for the critical exponents as a function of the control parameter a. a) The exponents /3 
(thin solid line), b = l3/u\\ (thin dashed line) and the mean-field prediction (thick solid line); note 
that (3^^ = 8^^^ . b) The dynamical exponent z = v^jv^ (thin solid line) and the mean-field predic- 
tion (thick solid line). The thick dashed lines on the right of both plots mark the best numerical 
estimates of the DP exponents in ID. 

C. Comparison with unrestricted Levy flights 



As anticipated in the Introduction, it is instructive to compare the results of the a-process 
with those of contact processes with infections mediated by unrestricted Levy flights (see 
Refs. jlfil. Ilfijl). In the latter case three different dynamical regimes were identified: 

• a mean-field regime for a < 3d/2 characterized by the critical exponents (3^^ = 1, 

= [a — d)~^, and = 1, where d denotes the spatial dimension. 

• a non-trivial phase for 3d/2 < a < a* with continuously varying exponents, restricted 
by the additional scaling relation z/y — ^±{0 — 2d) —2/5 = 0. 

• a DP regime for a > a*. 

Here the upper threshold for a can be expressed in terms of DP exponents by a* = 2d+z^^ — 
2(3°^ /v^ ~ 3.0766(2). The phase structure of this model exhibits a mean field regime for 
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a < 3d/2 which is not present in the a-process. In this mean field regime the unrestricted 
Levy fiights become so long-ranged that spatial fiuctuations are destroyed, leading to a 
homogeneous distribution of active sites. In the present model, however, the long-range 
interactions are effectively cut off at the distance to the nearest active site, hence there is 
no direct communication by long-range interactions between adjacent intervals of inactive 
sites. This restriction yields non-trivial correlations. 

It is instructive to compare mean-field equation with the corresponding equation for 
unrestricted Levy fiights, which - according to the present notation - can be written as 

dtp{^, t) = np{x, t) - p2(x, t) + DV°-V(x, t). (16) 

Here V"~^ is a symetric fractional derivative, which generates unrestricted Levy fiights. 
Clearly, this equation has scaling properties different from Q . For instance, within the 
mean-field approximation the exponent /3 is equal to 1. Therefore, it is likely(?) that also 
in low dimensions models with truncated and unrestricted Levy fiights belong to different 
universality classes. 

D. The a-process in higher dimensions 

The a-process and its mean-field equation can be easily generalized to d > 1. In this case 
the range, /, of the activation process of an inactive site located at i is given by the radius of 
the maximal "empty" (inactive) sphere centered in i. Within the mean-field approximation 
the average volume of the "empty" spheres is proportional to 1/p, so that the average radius 
of empty spheres is < / >oc p^^^'^. Moreover, the diffusion coefficient -D(p) can be easily 
computed by extending Eq. (jHl) to arbitrary dimension d 

D{p) ~ p^^-i + Do (17) 

Therefore, the corresponding mean-field equation reads 

dtp = rp- up'^l'' - V + ^(P) V'p + o(p') , (18) 

where r(g) vanishes at the critical point while M(g) and f (g) are positive constants. As in 
the \d case, \i a < d the absorbing state cannot be reached and there is no transition. On 
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the other hand, for a > 2d and d >2 a DP regime sets in, where both leading terms in the 
dynamical equation and in the diffusion constant become DP-like. In the range d < a < 2d 
the system displays a non-trivial transition, with continuously varying critical exponents. 
In summary, the mean-field exponents for d > 2 are 

z/f = 1 , (19) 

I ,\ -, d < a < 2d 

^MF ^ ^MF^MF ^ ^ a/d-1 ^20) 

1 a> 2d 



and 



^11 , d < a < d + 2 

MF II J a—d 

i a>d+2 



(21) 



III. GENERALIZATION OF THE MODEL 
A. The (T-process 

In a recent paper we studied a different version of a contact process with long-range 
interactions, called a-process 2^. In this model infections are short-ranged, i.e., active sites 
can only activate their nearest neighbors, but the rate for short-range infections depends 
algebraically on the distance m to the nearest active site before the update. Specifically, the 
cr-process is defined by the transition rates 

1^0 with rate 1 (22) 
10,01 11 with rate g(l + a/m'^) (23) 

where m denotes the distance to the next active site. Here a is a constant, q is again the 
control parameter, and a is an exponent controlling the characteristics of the interaction. 
It turns out that for a > 1 the transition belongs to the DP universality class, while for 
< o" < 1 the transition becomes first order. 

The model was introduced as a toy model in order to explain why non-equilibrium wet- 
ting processes, with a sufficiently strong attractive short-range force between interface and 
substrate, can exhibit a first-order transition. In this case the growth rate in the bulk is 
positive, so that spontaneous pinning far away from the edges can be neglected, meaning 
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that the transition is driven by unpinning and spontaneous growth of bound regions. The 
binding rate, however, was found to vary with the actual size of the unpinned regions, which 
motivates the algebraic m-dependent rate in Eq. (j2SI)- It should be noted that the cr-process 
is mainly relevant for non-equilibrium wetting in one spatial dimension only, where pinned 
sites effectively separate depinned islands into non communicating regions. Moreover, there 
is no unique and straightforward way to extend such a process in higher dimensions. In 
the following we combine the two long-range processes in a single Id model and study the 
corresponding phase diagram in the a-cr-plane. 



B. Definition of the a-a-process 

The combined a — a model, which evolves by random sequential update over a time step 
dt = 1/L is defined by the following transition rates: 

1^0 with rate 1 (24) 
10,01 11 with rate g(l + a/m'") (25) 
000 ^ 010 with rate bq/T (26) 

where m is the size of the inactive island and / denotes the distance from the nearest active 
site. Here q is the control parameter, a and a are the control exponents, and a and b are 
constants. The "pure" a- and a- processes are recovered in the limits a ^ oo and a oo, 
respectively. 



C. Mean-field analysis and numerical simulations 



Neglecting spatial 



(see Sec. Ill Al and |^|) becomes 



fluctuations, the mean field equation of the combined a — a-process 



d,p = -p + qp' + pT + &?P(2 -P)Y. ^ 'I^J . (27) 

m=l 1=2 

In the regime of interest, namely for < a < 1 and 1 < a < 2, the mean-field equation can 
be written to leading order as 

dtp = rp+ pp'+'' - up'' + 0(p2) , (28) 
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where r, p and u are g-dependent constants: 



u 



P 



r 



q — 1 + bq 



bq -r(2-a). 

a — i 




a 



- 1 



(29) 
(30) 
(31) 



Note that these are approximate expressions, as they are obtained by substituting the 
sums in ()27|) by integrals to obtain simple analytic expressions. This approximation does 
not affect the sign of p and u which turn out to be positive constants in the regime of 
interest. Accordingly, if the leading nonlinear term is p^^'^ (i.e. a < a — 1) , the transition 
is first order. On the other hand, if the leading nonlinear term is p° (cr > a — 1) , the 
critical behavior of the pure a-model is recovered. Thus, the mean-field approach predicts 
that the critical exponents are independent of cr . One can add the effect of the long range 
processes on the diffusion constant as was done in Sec. Ill Al It is readily seen that the 
critical exponent z assumes its DP value for a > 3 as long as a > 1. 

The mean-field phase diagram is shown in Fig. |3] It contains four different phases: 

(i) a DP phase for cr > 1, a > 3, 

(ii) a continuously varying exponent phase CVEl for a > a — 1 and 1 < a < 2 where the 
transition is second order with continuously varying exponents (3, 6 and z. 

(iii) a continuously varying exponent phase CVE2 for cr > 1 and 2 < a < 3 where the 
transition is second order but where z is the only continuously varying exponent. 

(iv) a phase for cr < a — 1 and < cr < 1 
where the transition is discontinuous. 

Numerical simulations have been performed for testing the predictions of the mean-field 
phase diagram. In these simulations we have taken a = 2 and 6=1. 

An extensive numerical study of the phase-diagram of the a-cr-process is time consuming 
and was not tackled in the present paper. Rather, we restricted our analysis to a small 
number of cuts in the a-a plane in order to test the validity of the mean-field predictions. 
We find that while the general features of the phase-diagram are correctly reproduced by 
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a 


2.0 


2.3 


2.5 


2.6 


2.7 


3.0 


(Ic 


0.7783(2) 


0.9308(3) 


1.0155(3) 


1.0530(3) 


1.0880(5) 


1.1785(5) 


Z 


1.18(3) 


1.40(3) 


1.51(3) 


1.50(4) 


1.54(5) 


1.57(3) 


5 


0.32(2) 


0.22(1) 


0.19(1) 


1.18(1) 


1.18(1) 


1.16(1) 


(3 


0.52(1) 


0.37(3) 


0.32(3) 


0.30(2) 


0.29(1) 


0.28(2) 



TABLE II: a — a process: Estimates of the critical points and exponents for a = 2 and various 
values of a. Other parameters have been fixed to a = 2 and 6 = 1 

the mean-field, the location of the transition lines and the values of the critical exponents 
are changed. 

The numerical method proposed in [2^ has been applied for the inspection of the phase 
diagram along the line a = 2 for a discrete set of values of a ranging between 2 and 3 . In 
these simulations we have used the same system sizes and conditions described in |n] The 
numerical simulations suggest that in fact there is a single phase with continuously varying 
exponents (CVE) rather than two as predicted by the mean field. In this phase all critical 
exponents are continuously varying with a. Both in the DP and CVE regions the critical 
exponents coincide, within numerical accuracy, with those found for the pure a-model (i.e., 
a ^ 00). The results reported in Table ITTl are in agreement with the mean-field prediction 
that scaling properties in the critical regions are independent of the a-process. It is worth 
stressing that even for finite values of cr the boundary between the DP and the CVE region 
is located close to a ~ 2.7 . The position of the boundary line for o" = 00 and cr = 2 is 
marked by crosses on the phase diagram (see Fig. ^ 

We now consider the boundary between the first-order and the DP regions. For the 
pure cr-model (i.e., a 00) this boundary was investigated numerically in Ref. [2^ . 
Following ( 2^) we determine the nature of the phase transition for finite a by studying 
the size distribution of inactive domains in stationary active state. At the transition to the 
absorbing state this distribution is expected to exhibit a power law tail P(m) ~ m"'*' for 
large inactive domains of size m. By numerically determining the exponent 7, the nature of 
the transition may be deduced. For 7 > 2 the average domain size is finite at the transition 
and thus it is first order. On the other hand for 7 < 2 the average domain size diverges 
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1.0184(1) 


1.0430(5) 


1.0647(3) 


1.0992(1) 


1.1445(5) 



TABLE III: a — cr-process: estimates of the critical points qc for Q! = 3, a = 2, & =1 and various 
values of o. 

at the transition and thus the transition is continuous. In particular, for DP, where the 
transition is continuous, one has 7°'' = 2 — /?/z/_l ~ 1.747 . 

We have analyzed the nature of the phase transition for a = 3 at some values of a ranging 
between 0.5 and 2 . We find that 7 = 7^^ for a > 0.8 indicating that the transition to the 
absorbing state is of DP nature. On the other hand for cr < 0.8 we find 7 > 2, indicating a 
first order transition. The position of the boundary line for a = cxo and a = 3 is marked by 
dots on the phase diagram (see Fig. H}. 

Details of the numerical analysis is presented in Fig. |3| In this figure we first illustrate the 
method for identifying the transition point. This is done by plotting p{t)t^^'' as a function 
of t for different values of q and searching for the q value for which this quantity approaches 
a constant. The data were obtained for for quite large system sizes (L ^ 2^*^ ~ 2^^) and 
homogeneous initial conditions. We then display the distribution function P{m) at the 
transition, and determine the exponent 7 which controls its large m behavior. The value of 
the control parameter q at the transition point is given in Table lllll for a = 3 and several 
values of the parameter a. 

Our numerical results suggest that the boundary between first order and DP like behavior 
varies with a. In particular it takes place at cr ~ 0.8 for a = 3 while it is a ^ 1 in the 
limit of infinite a. Thus the boundary line in the a — a plane is not parallel to the a axis 
at variance with mean field prediction. 

IV. CONCLUSIONS 

In many cases it is useful to project complex dynamical processes onto simpler ones which 
could be more readily analyzed. In doing so it may happen that the local dynamics of the 




FIG. 4: Mean-field phase diagram of the a-a-process. Four different regions are identified corre- 
sponding to different types of transitions to the absorbing state: a DP region, a region of a first 
order transition and two regions of continuously varying critical exponents as explained in the text. 
We also indicate the of the phase boundary lines as determined by numerical simulations: full dots 
for the DP-first order boundary and crosses for the DP-CVE boundary. The numerical results 
suggest the existence of a single CVE phase with a single DP-CVE line located at a ~ 2.7 . 



original model is translated into a non-local dynamics of the projected model. For example 
it has been argued that when the dynamics of fluctuating interfaces interacting with a wall 
is projected onto the dynamics of a contact process different types of long range processes 
may take place. 

Motivated by this observation we introduced a non-local version of the contact process 
where particles are created with a rate which decays like with the distance / from the 
nearest particle. Mean-field analysis shows that the critical exponent (3 associated with the 
transition to the absorbing state varies continuously with a for 1 < a < 2 and it diverges 
for a approaching 1. The usual short range mean-field exponent /? = 1 is recovered above 
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FIG. 5: Numerical study of the a-cj-process for a = 3.0. a) p{t)t vs t for a = 
0.8 and values of the control parameter close to qc ~ 1.043 : from top to bottom q = 
1.045, 1.044, 1.043, 1.042, 1040 . b) stationary density pstat as a function of g for a = 0.7: 
a jump occurs at the first-order transition point qc ~ 1.0184 . In c) and d) we display the proba- 
bility distribution of the size of inactive domains P{m) at the transition, c) from top to bottom 
a = 1.5, 1.1, 0.9, 0.8) : the power-law decay agrees with the DP scaling (see dashed line), d) from 
top to bottom (T = 0.7, 0.5 : P{m) decays faster then so that (m) is finite. In both of these 
plots the curves are not normalized and have been rescaled to fit on the same plot. 

the threshold a = 2. The exponent v\\ is found to be 1 in the entire range. Numerical 
studies in one dimension support this general behavior, although the value of (3 and the 
upper threshold seem to be different. 

By including spatial fluctuations in the mean-field equations the mean-field exponents v± 
and z are evaluated. It is found that near a = 1 the exponent z vanishes linearly, while it 
approaches the usual mean-field value z = 2 at another threshold, which in one dimension 
is a = 3. Also in this case we find a qualitative agreement with numerics. 
Extension of the model to include another long range dynamical process, the cr-process, has 
also been considered using both mean-field approximation and numerical simulations. Here, 
too, the global features of the mean-field phase diagram are qualitatively recovered by the 
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numerical study. 
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